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We have studied magnetic vortex oscillations in soft sub-micron cylindrical dots with variable 
thickness and diameter by an analytical approach and micromagnetic simulations. We have con- 
sidered two kinds of modes of the vortex magnetization oscillations: 1) low-frequency translation 
mode, corresponding to the movement of the vortex as a whole near its equilibrium position; 2) high- 
frequency vortex modes, which correspond to radially symmetric oscillations of the vortex magnetiza- 
tion, mainly outside the vortex core. The vortex translational eigenmode was calculated numerically 
in frequency and time domains for different dot aspect ratios. To describe the discrete set of vortex 
high-frequency modes we applied the linearized equation of motion of dynamic magnetization over 
the vortex ground state. We considered only radially symmetric magnetization oscillations modes. 
The eigenfrequencies of both kinds of excitation modes are determined by magnetostatic interac- 
tions. They are proportional to the thickness/diameter ratio and lie in the GHz range for typical 
dot sizes. 

PACS numbers: 75.40.Gb,75.30.Ds,75.50.Bb,75.75.+a 



I. INTRODUCTION 

The magnetization vector of a ferromagnetic body 
(particle) oscillates when variable external magnetic 
fields are applied. The amplitude of magnetization os- 
cillations, which depends on the intrinsic fields as well as 
on the particle's geometrical parameters increases dras- 
tically when the frequency of the external field coin- 
cides with the frequency of a fundamental spin-excitation 
eigenmode of the system. 

This issue has important implications for developing 
high-speed and high-density magnetic devices, where 
short field pulses are applied to magnetic media to 
achieve the magnetization reversal— The characteristic 
switching time is usually in the nanosecond and sub- 
nanosecond range for non-uniform and coherent magne- 
tization rotation, respectively. Understanding of spin- 
wave excitation spectra is thereby crucial to determine 
the field-dependent spin-instability regions, where spon- 
taneous or thermally-assisted magnetization reversal pro- 
cesses might occur—. 

In the general case, both the short-range exchange and 
the long-range magnetostatic interactions contribute to 
eigenfrequencies of the collective spin excitations. The 
long-wavelength excitations are mainly influenced by 
dipole-dipole interactions. Therefore, the demagnetizing 
fields of the magnetic elements determine their eigenfre- 
quencies. These fields can be strongly non-uniform in 
the case of non-ellipsoidal samples or if the static magne- 
tization distribution differs significantly from the single- 
domain (collinear) state. This is especially the case for 
technologically important flat magnetic particles, pre- 
pared by means of micro-fabrication and thin film de- 
position techniques. The adequate description of high- 
frequency magnetization dynamics in such systems is a 
challenge for modern magnetism theory. 

The remanent magnetization distribution within a fer- 



romagnetic particle is a complex phenomenon that de- 
pends on the size and shape of the particle as well 
as on the balance of different contributions (exchange, 
magnetostatic, anisotropy etc.) to the magnetic energy. 
Exchange energy dominates at small particle sizes and 
favors a uniform magnetization distribution, so-called 
single-domain or "flower" ("leaf") state. On the other 
hand, this almost collinear spin alignment leads to large 
demagnetizing fields due to magnetic surface charges and 
correspondingly increase of the magnetostatic energy. 
The importance of the magnetostatic energy increases 
gradually as the particle size increases. The competition 
between these two energies and anisotropy energy leads 
to non-uniform magnetic states such as multi-domain 
states (if the anisotropy constant is greater than the mag- 
netostatic energy) or magnetic vortices (if the anisotropy 
energy is small in comparison with the magnetostatic en- 
ergy) . The critical size at which a particle becomes single 
domain depends on the interplay between the energies 
mentioned above, and also on the shape of the particle. 
This critical size is of the order of the material exchange 
length L ox = y / 2A/M^ , where A is the exchange stiff- 
ness constant and M s is the saturation magnetization. A 
magnetically soft particle with lateral dimensions smaller 
than the exchange length (L cx fss 10 — 20 nm) is expected 
to remain in a single domain state at remanencc. 

Sub-micron magnetic elements with regular shapes 
(rectangular, cylindrical, spherical etc.) have attracted 
much attention during the last few years due to consid- 
erable progress in fabrication and characterization of the 
samples 

In this paper we will consider dynamic properties of 
circular particles (dots) made of soft magnetic material 
with a thickness of about L ex and diameters considerably 
larger than L cx . It is well established, that a magnetic 
vortex structure is the ground state of such submicron 
d t Sl 5j6 J 7 i 8j9 i io i ii This vortex state can be interpreted as a 
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part of a two-dimensional magnetic topological soliton 12 
adapted to the particle shape to minimize the total mag- 
netic energy. Experimental studies of the vortex state in 
flat submicron dots have been performed using different 
methods. Magnetization reversal in sub-micron dots due 
to magnetic vortex formation and its displacement has 
been explored by analytical and numerical micromag- 
netic modeling^ii^ It was found from static measure- 
ments and calculations that the magnetic vortex state 
corresponds to a deep energy minimum for submicron 
dot radii and it can be destabilized either by reducing 
the dot radius down to the material's exchange length 
or by applying an in-plane magnetic field. The vortex 
state is not specific for the circular shape, it was also ob- 
served in rectangular^ ellipticalji£*i£ and triangular 
flat magnetic elements. 

The dynamic properties of sub-micron dots have only 
more recently attracted some attention. 10 There have 
been a few successful experimental investigations of the 
magnetization dynamics under short field pulses in a sat- 
urated FeNi diskiS and also closure-domains in a Co 
disk^S In this article we concentrate on the linear spin 
dynamics in an equilibrium, but essentially non-uniform, 
vortex state. Spin excitations in the vortex-state are ex- 
pected to be substantially different from those in the uni- 
formly magnetized state and also from the excitations 
(spin-waves) observed in infinite thin magnetic films. 
An example is the appearance of a low-frequency mode 
(translation mode) associated with the movement of the 
vortex as a whole in a potential well created mainly by 
magnetostatic energyi^LSSi It was found that the vortex 
core undergoes a spiral motion with the eigenfrequcncy 
in the sub-GHz range for thin sub-micron dots. This 
lowest frequency translation mode can be treated as vor- 
tex core excitation and is localized near the dot center. 
This mode is similar to crossed-domain wall resonance 
in a magnetic film. 23 We consider that the spectrum of 
vortex spin excitations will be quantized due to parti- 
cle's (dot's) restricted geometry as was shown in 2 ^ 2 ^ for 
the saturated state. Spin excitations of the vortex state 
have been discussed for an infinite film (Ref. l26l and Refs. 
therein) neglecting the magnetostatic interaction. This 
is an unsuitable approximation to describe magnetization 
dynamics of sub-micron magnetic particles, where the ex- 
change contribution to the low-lying eigenfrequencies is 
essentially smaller than the magnetostatic contribution 
(L ex *C R, where R is the disk radius). Magnetization 
dynamics in this case is governed by the long-range dipo- 
lar forces. We are not aware of any theoretical or experi- 
mental investigation of dipole dominated high frequency 
spin excitations in a magnetic vortex state except recent 
results by Novosad et alii - and Ivanov et ali22 

In this paper we address the precessional spin dynam- 
ics of the vortex ground state and determine the low- 
lying vortex eigenfrequencies which correspond to radi- 
ally symmetric spin excitations. The eigenfrequencies are 
well above the vortex translation eigenfrequency. Analyt- 
ical calculations of the eigenfrequencies are supported by 



numerical micromagnetic simulations. 

The vortex magnetization dynamics can be detected, 
for instance, by ferromagnetic resonance (FMR), Bril- 
louin light scattering (BLS) techniques or by time- 
resolved Kerr microscopy. The spectrum of the dot dis- 
crete eigenfrequencies can be used as a starting point for 
the consideration of magnetization reversal. 

We present in Sec. |H] a theoretical description and in 
Sec. IIIII micromagnetic calculations of discrete eigenfre- 
quencies of excitation modes of the vortex state in soft 
magnetic cylindrical dots. We consider that above the 
vortex translation mode frequency (which is < 1 GHz) 
there is a set of radially and azimuthal symmetric modes 
(of about 10 GHz) localized outside the vortex core. For 
sub-micron dot sizes the dynamic demagnetizing fields 
determine the eigenfrequencies of these modes, whereas 
the translation mode frequency is determined by the 
static demagnetizing field (magnetostatic contribution to 
the restoring force) . The vortex excitation eigenfrequen- 
cies are found from eigenvalues of a magnetostatic inte- 
gral operator (Sec. [HJ and also from dynamical micro- 
magnetic simulations (Sec. IIIII) based on the Landau- 
Lifshitz equation of motion. Full details of the ana- 
lytical model are reported. Experimental techniques to 
study these vortex modes are discussed. A comparison 
of numerical (FEM-BEM and finite differences) micro- 
magnetic simulations with analytical results is presented. 
Finally, the summary is given in Sec. IIVI 



II. THEORY 

We present a theoretical description of the vortex state 
oscillations in submicron-size cylindrical dots. Our ap- 
proach is based on the consideration of small dynamic 
oscillations over the centered vortex ground state. We 
consider here only radially symmetric vortex excitation 
modes in zero applied magnetic field. These modes have 
frequencies well above the vortex core translation fre- 
quency for the given dot sizes. They have maximum 
FMR and BLS intensity due to their non-zero average 
dynamic magnetization. It is known that the vortex core 
bears some topological charges (gyrovector and vorticity, 
see Ref. Il2j) and this peculiarity has important conse- 
quences for the vortex state dynamics. In the follow- 
ing, we do not consider the vortex charges and assume 
that the high-frequency excitation modes are localized 
outside the vortex core. In other words, we neglect the 
vortex core influence on these modes. Let L denote the 
dot thickness, and R the dot radius. Since the dot thick- 
ness is small, we assume a two-dimensional magnetization 
distribution m = M/A/ s , m 2 = 1 within the cylindrical 
dot, which does not depend on the z-coordinate along 
the dot thickness. We use the angular parameteriza- 
tion for the dot magnetization components m z = cos0, 
m x + irriy = exp(i$)sinO, $ — </> ± n/2 . Here posi- 
tive and negative signs correspond to counter-clockwise 
and clockwise rotation of the vector m in the dot plane, 
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respectively. The cylindrical coordinates p, <f>, z are used 
below. The spin structure of the static vortex located in 
the dot center may be described by the one-parameter 
ansatz 



(JJL)V(/0= / dp'p'r(p,p')fj,( P '), «,) 



tam(e (p)/2) = p/b 



(1) 



if p < b and 80 = tt/2 if p > b as suggested by Usov 
et ali££ The radius of the vortex "core" b « 10 nm can 
be determined from magnetic energy minimization. The 
total dot magnetic energy including exchange and mag- 
netostatic energies is: 



W = L d 2 p[A(VQ) 2 + sin 2 6(V$) 2 ) - -H m ■ M] 



(2) 

where H m is the magnetostatic field. We use the equa- 
tion of motion of variable magnetization 
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The calculation of the exchange contribution to the 
effective field is straightforward, while the magneto- 
static field is expressed via the magnetostatic Green's 
function 24 in cylindrical coordinates p,(j),z. The mag- 
netization distribution m(r) within the cylindrical dot 
leads to the variable demagnetizing field h(r) = G[m(r)], 
where G is a tensorial non-local integral operator (the 
tensorial magnetostatic Green's function G a p(r, r')) ex- 
pressed in the cylindrical coordinates 



G[m(r)] = / G[r,r']m(r')dV, 
G Q/3 (r,r') = -(V r )a(V P /) / s 1 ' 



r — r 



(4) 



We use averaging over z, z' and also over 4>' (due to 
the assumed radial symmetry of m(r)). Only pp and zz- 
components of the averaged G tensor are not equal to 
zero. The equation of motion |J3J| is then linearized by 
substituting m(p,t) — m (p) + p(p,t) or by magnetiza- 
tion angles: 



G(p, t) = 9 (p)+i%, t), $(p, <p, t) = v +-+^(p, t), (5) 

where the vector p = {pp,p Vl p z ) = 
(— ^ sin 0q, I? cos 0o, — $sin0o) is the dynamic part 
of the magnetization. For small variables d and 
p = p p = —ip sin Qq we get a system of two coupled 
integro-differcntial equations. We can neglect the 
dynamic exchange interaction for sub-micron dot radii, 
exclude the variable $(p, t) and reduce the problem to 
an eigenvalue problem for the integral magnetostatic 
operator: 



T (P,P')= / drrsm@o(r)gzz(p,r)g P p(r,p'), (7) 



where the averaged over thickness Green's function com- 
ponents do not depend on 



9afl{p,P ) = ~j- 
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Solution of Eq. © yields the cigenfunctions p n (p) and 
corresponding eigenfrequencies oj n , which can be num- 
bered by integers n — 1,2,... omitting the azimuthal in- 
dex m = (radial symmetry). We use the index n = for 
the vortex translation modef^i This mode has azimuthal 
indices m = ±1. We do not assume any boundary con- 
ditions for the dynamic magnetization on the dot side 
boundary due to negligibly small surface anisotropy for 
soft magnetic materials (FeNi, for instance). The vortex 
core is a stable formation which is strongly coupled by 
exchange forces and we assume that the dynamic magne- 
tization has considerable values only outside of the vor- 
tex core, in the area where &o(p) = tt/2 . Then the 
dynamic magnetization will have only two components, 
p p and p z given by Eq. (JJJ, which are perpendicular to 
the vortex static magnetization mo = (0, mg,0). These 
components depend on time as p p (t) ~ sin(wt + 4>vi), 
P-z(t) ~ cos(ujt + <po), where </>o is an initial phase. This 
time dependence corresponds to elliptical precession of 
the dynamic magnetization components in the p — z 
plane around the vortex static equilibrium magnetiza- 
tion mo. The dynamic p p component determines the 
dynamic side surface and volume magnetic charges. The 
coupled oscillations of p p and p z then lead to a con- 
siderable increase of the dynamic magnetostatic energy 
(mainly due to volume charges, because surface charges 
described by p p (R) magnetization component are small) 
and to more high eigenfrequencies in comparison with 
the case of the vortex translation mode (no induced sur- 
face charges). Assuming that the dynamic components 
are concentrated mainly outside of the vortex core, in the 
area &o(p) = 7r/2, for small dot aspect ratio L/R <C 1 
(in this case g Z z(p,p') — — ^5(p — p')/p) we get from 
Eqs. and |JSJ the simplified integral equation 
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The eigenfrequency w„ which corresponds to the normal- 
ized to unit n-eigenfunction p n (p) is 



CD 
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47T 



dpp / dp' p' g pp {p, p')p n {p)p n {p') 
Jo 

(11) 

In this magnetostatic approximation the vortex eigen- 
frequencies depend only on the dot aspect ratio L/R. 
More careful analysis shows that the eigenfrequencies are 
approximately proportional to \J L/R if the aspect ratio 
(3 = L/R -C 1. The solution of the spectral problem 
Eq. ® will give us a discrete set of approximate magne- 
tostatic eigcnfunctions (radial modes profiles) and eigen- 
frequencies. The kernel of Eq. © is real and can be 
symmetrized. This means that the Hilbert-Schmidt the- 
ory of the integral equations is applicable for our case. 
The eigenvalues are real and the corresponding eigen- 
functions are orthogonal. The approximate numerical 
solution of the spectral problem (Eq. (0) for FeNi dot 
parameters (with M s = 800 G, j/2n = 2.95 GHz/kOe 
and wm/27t = 29.7 GHz) and different L/R gives us the 
set of eigenfrequencies and eigenfunctions (see Fig. ^ HI 
and 0J) ■ In the limiting case of a thin dot j3 <C 1 this 
problem can be solved analytically. We get eigenfrequen- 
cies (w„/wm) 2 = f{fiot n ) and eigenfunctions p n {p) — 
C n Ji(K n p), where a n is the nth root of the equation 
J\{x) = 0, K n = a n /R, Ji(x) is the first order Bessel 
function, and C n is the normalization constant. The set 
of be treated as a set of allowed quantized radial 

wave numbers, which is consistent with the geometry of a 
circular dot and its symmetry of the vortex ground state. 
The equation for n n corresponds to the strong pinning on 
the dot side surface p = R. In contrast, the numerical so- 
lution of Eq. 10 corresponds to some intermediate values 
of the pinning, although this pinning is strong enough. 
This pinning is of pure dipolar origin and it was calcu- 
lated for thin magnetic stripes in Ref. |23- The radial 
eigenfunctions Ji(n n p) form a complete orthogonal set 
and can be considered as good trial functions, which are 
asymptotically exact within the limit — » 0. The eigen- 
frequencies u) n correspond to dynamic magnetostatic os- 
cillations in the absence of the internal magnetic field. 



III. MICROMAGNETIC CALCULATIONS AND 
DISCUSSION 

A. Vortex static configuration 

The numerical 3D dynamic micromagnetic simulations 
are carried out using a hybrid finite element /boundary el- 
ement method.— The nanodots are discretized into tetra- 
hedral finite elements. The magnetic polarization is de- 
fined on the nodes of the finite element mesh and lin- 
ear test functions are used on isoparametric elements. 
The effective field which is the sum of the exchange, 
anisotropy, magnetostatic, and external field is calcu- 
lated directly (for the local contributions from exchange, 
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FIG. 1: Dependence of the eigenfrequencies of the vortex ra- 
dial modes with the indices n=l,2,3 on the dot aspect ratio 
L/R. The closed symbols correspond to numerical solutions 
of Eqs. and itTTJt. 
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FIG. 2: The vortex radial eigenfunctions Hn(p) for the dot 
aspect ratio (3 = 0.075 calculated by Eqs. © and 1101 . The 
solid lines are numerical solutions of Eqs. @ and 1101 . the 
dashed lines are the Bessel functions Ji(k„p). 
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FIG. 3: The vortex radial eigenfunctions p n (p) for the dot 
aspect ratio f3 = 0.15. The solid lines are numerical solutions 
of Eqs. © and 1101 . the dashed lines are the Bessel functions 
Ji(K n p). 
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FIG. 4: The calculated profile of the out-of-plane M z mag- 
netization component in a circular nanodot (R = 200 nm, 
L — 10 nm). The solid line is solution of Eq. the dashed 
line corresponds to finite element micromagnetic simulations. 

anisotropy, and external field) as H c ff ~ —dW/dm (W 
is the energy density in Eq. and using a boundary 
element method^ for the long range magnetostatic field 
with open boundary conditions. The dynamic Landau- 
Lifshitz equation of motion for the magnetization is in- 
tegrated with a preconditioned backward differentiation 
methodi^Li^ For the numerical micromagnetic finite el- 
ement simulations we have assumed the material pa- 
rameters of permalloy (NisoFe2o), which is a typical 
soft magnetic material with negligible magnetocrystallinc 
anisotropy. We set the magnetic saturation polarization 
to J s = fJ,oM s « 1 T and the exchange constant to 
A = 13 x 10~ 12 J/m, which gives an exchange length 
L ox = y / 2A/M2 = 5.7 nm. 

The finite element simulations have been initialized 
with the magnetization distribution of the vortex model 
based on Eq. (JJ) and an approximate core radius. Then 
the Landau-Lifshitz equation of motion for the magne- 
tization has been integrated with a damping constant 
a = 1 in zero field, and the magnetization relaxed to 
its equilibrium distribution, which minimizes the total 
Gibbs' free energy. The plots of the normal magnetiza- 
tion component M z in Fig. 0] show a comparison of the 
analytic vortex model with the equilibrium magnetiza- 
tion distribution of the FE simulation. 

We define the vortex core radius by the equation 
Mz{p) = 0. The numerical results show, that the vortex 
core is approximately 50% larger (18.5 nm) than assumed 
by the vortex model given by Eq. Q (12 nm). Further- 
more the finite element simulation shows that there is 
a region with M z < outside the core. Thus, we find 
positive surface charges in the core of the vortex, which 
are surrounded by negative surface charges. Only outside 
of approximately half the radius (50 nm) almost all sur- 
face charges disappear. It has been verified, that there 
is very little variation of the magnetization distribution 
across the thickness of the nanodot. The effect of small 
oscillations of the static vortex magnetization profile can 
be neglected for large R > 200 nm. 

The difference in magnetostatic and exchange energy 
between the approximate analytical vortex model given 
by Eq. and the finite element simulation is calculated. 
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FIG. 5: Oscillations of the in-plane averaged vortex mag- 
netization component (M x ) as a function of time (a) and its 
Fourier spectrum (b) for the dot with L/R — 20 nm/100 nm = 
0.2 under applied in-plane field H x — 0.01 T. 



The static vortex magnetization distribution within the 
analytic vortex model has been obtained using the formu- 
las given in Refs.l27landl33L The equilibrium magnetiza- 
tion distribution after relaxation (dashed line in Fig. 
as calculated with the finite element model leads to a 
4.4% lower total energy. We confirmed numerically that 
the vortex core has relatively small radius < 20 nm and, 
therefore, it can be ignored in the description of magne- 
tization dynamics of large dots with 0.2 < R < 2/j,m£2^ 



B. Vortex translation mode 

To study the vortex lowest frequency mode (transla- 
tion mode) numerically an external field of 10 mT has 
been applied in the plane of the dot to the remanent 
state and a damping constant of a = 0.05 has been as- 
sumed. The time evolution of (M x ) (the average of M x 
over the whole nanodot) for a dot with an aspect ratio of 
L/R = 10 nm/100 nm = 0.1 is given in Fig.GJa). Then 
the damped oscillation, which is caused by the spiral mo- 
tion of the vortex core towards its equilibrium position, is 
observed. The corresponding Fourier spectrum is given 
in Fig. E£b) and shows a sharp peak at a frequency of 
0.7 GHz. 

Fig. [f)| shows the results of calculations of the trans- 
lation mode eigenfrequencies of various nanodots with 
a radius R = 100 nm and a thickness between 10 and 
40 nm. The translation eigenfrequency is approximately 
linearly proportional to the dot aspect ratio L/R. The 
results are in good agreement with the results of a finite 
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FIG. 6: Translation mode eigenfrequency versus the dot as- 
pect ratio L/R for cylindrical nanodots with 7? = 100 nm. 
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FIG. 7: The contributions to magnetic energy vs. time for 
a cylindrical dot with L/R = 0.1 and H x = 0.01 T. Solid 
line: the total magnetic energy, dashed line: exchange energy, 
dotted line: magnetostatic energy, dot-dashed line: Zeeman 
energy. 



difference model and the analytical "two-vortices" model 
presented in Ref. |2lJ 

The decreasing total energy (dissipation due to damp- 
ing with a = 0.05 in the Landau-Lifshitz equation of 
motion) and the swapping between magnetostatic and 
Zeeman energy are shown in Fig. [7] The exchange en- 
ergy remains constant, because the vortex core, which 
accounts for most of the exchange energy, moves without 
changing its shape. This supports the analytical descrip- 
tion of the translational mode suggested in Ref. 21, where 
circular motion of the vortex core as whole was calcu- 
lated. Note that recently existence of this low-frequency 
vortex translation mode was confirmed experimentally 
by direct Kerr microscopy measurements in Ref. l35l The 
eigenfrequency measured by Park et al. agrees with our 
calculations with an error of less than 20%. 



C. Vortex radial modes 

The vortex radial modes were excited by applying 
an external field of 5 mT in z-direction and relaxing 
the magnetization with large damping (a = 1). When 
equilibrium was reached, the external field was instanta- 
neously switched off and the free oscillation in zero field 
with very low damping was studied. With a = 0.05 the 
oscillations are damped out within a few oscillations. As 
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FIG. 8: Oscillations of the averaged out-of-plane magnetiza- 
tion component (M z ) as a function of time (a) and its Fourier 
spectrum (b) for the cylindrical dot with L/R = 0.2 under 
applied out-of-plane held H z — 0.005 T. 



a result the Fourier spectra are very poor because the 
resolution of the Fourier spectra increases with measure- 
ment time. Therefore a = 0.001 has been used to observe 
many oscillations and obtain a high resolution Fourier 
spectrum. 

The time dependence of the average magnetization 
(M z ) and the Fourier spectra for nanodots with an as- 
pect ratio of L/R = 20 nm/100 nm = 0.2 and L/R = 
40 nm/200 nm = 0.2 are given in Figs. Eta) and (b) 
and Eta) and (b), respectively. Note that the azimuthal 
modes (with the index to ^ 0) do not contribute to oscil- 
lations of (M z ). We checked numerically that the calcu- 
lated modes have radial symmetry and correspond to the 
rotation of the dynamical magnetization [i = (/z p ,0,/z z ) 
around the vortex static magnetization mo = (0, TOq,0) 
(cf. Fig. EE). 

For the constant aspect ratio L/R = 0.2 we find an 
eigenfrequency of approximately 12.6 GHz. However, for 
a nanodot with L/R = 40 nm/200 nm = 0.2 a very pro- 
nounced beating is observed (Fig. |SJ. This is due to the 
fact, that there is another eigenfrequency of 11.6 GHz 
very close to the 12.6 GHz oscillation. However, the 
main peak position depends only on the combination ra- 
tio L/R. This confirms the magnetostatic origin of the 
mode. But the physical picture is more complicated for 
large L, when the magnetization dependence on z may 
be essential. 

We concentrated on numerical simulations of the de- 
pendence of the first magnetostatic radial eigenfrequency 
on the dot aspect ratio L/R. The calculated evolution 
of the first radial mode with increasing dot thickness L 
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FIG. 9: Oscillations of the averaged out-of-plane magne- 
tization component (M z ) as a function of time (a) and its 
Fourier spectrum (b) for the cylindrical dot with L/R — 
40 nm/200 nm = 0.2. 
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FIG. 10: Oscillations of the vortex magnetization compo- 
nents, (a) M z (measured at different positions £ = x/R, solid 
line: £ = —1, dotted line: £ = —0.8, dashed line: £ = —0.5, 
dot-dashed line: £ = +0.5) oscillates in phase across the whole 
nanodot. The oscillation amplitude slightly decreases towards 
the center, (b) M p and M z as a function of time at different 
radii r = p/R. Solid line: M z at r — 1.0, dotted line: M p 
at r = 1.0, dashed line: M z at r = 0.5, dot-dashed line: M p 
at r = 0.5). The phase shift between the M p and M z com- 
ponents is equal to 90 deg (magnetization rotation in the z-p 
plane) . 
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FIG. 11: Evolution of the Fourier spectrum of the component 
(rn z ) for the first (n = 1) radial vortex excitation mode with 
varying the dot thickness L from 5 nm to 30 nm. R — 200 nm. 
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FIG. 12: The radial vortex mode eigenfrequencies versus as- 
pect ratio L/R for nanodots with different thickness L. The 
open triangles are points simulated numerically, the solid line 
corresponds to fitting using the equation cv n ~ y L/R. 



(5-30 nm) is shown in Fig. 1111 The eigenfrequency of this 
mode increases with L in agreement with analytical cal- 
culations of the vortex magnetostatic modes in Sec. [H] 
We plotted the micromagnetically simulated eigenfre- 
quency as a function of the dot aspect ratio L/Rin Figll2l 
together with the analytic dependence u> n ~ y/ L/R and 
got a good agreement of the eigenfrequencies calculated 
by these different methods. Thus, our numerical cal- 
culations confirmed that the radial eigenfunctions and 
eigenfrequencies which correspond to low-lying part of 
the vortex dot excitation spectrum have magnetostatic 
origin. 

The calculated radial eigenmodes are magnetostatic 
(wave length is of the order of the dot radius) and concen- 
trated in the area outside of the vortex core (see Figs.|2E] 
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for the calculated eigenfunctions). The eigenfrequencies 
of these modes u> n are determined by the dynamic de- 
magnetizing fields. The static dipolar field for the vortex 
structure has two components H z and H p , which depend 
on p only due to circular symmetry of the vortex. Both 
components are negligibly small in the area of our inter- 
est, i.e. outside the vortex core. In this sense the problem 
of non-uniform magnetostatic spin-excitations is similar 
to the problem of magnetostatic waves in thin magnetic 
stripes magnetized alon g th eir long side considered by 
Guslienko et al. in Ref. l28l Effective boundary condi- 
tions for dynamic magnetization on the dot surface side 
can be also formulated similarly to Ref. |28j The condi- 
tion Hn{p = 0) = is satisfied near the dot center. This 
is consistent with our assumption about concentration 
of the dynamic magnetization of the radial modes out- 
side the vortex core. These radially symmetric magnetic 
dynamic excitations can be detected by modern experi- 
mental techniques such as BLS, FMR,and time-resolved 
Kerr- and magnetic circular dichroism measurements. 

Very recently, when this manuscript was ready for sub- 
mission, we became aware of experimental data by Buess 
et ali^ii They measured the eigenfrequencies of the vor- 
tex radial modes by time-resolved Kerr microscopy and 
obtained frequencies of 2.8, 3.9 and 4.5 GHz. The dot 
aspect ratio (3 = L/R = 15 nm/3000 nm = 0.005 is 
very small and we can apply the approximate equation 
(uj n /idM) 2 = f{fia n ) derived from Eq. (JTTJ. The calcu- 
lated eigenfrequencies 2.91, 3.93 and 4.73 GHz are very 
close to the experimental data. 



IV. CONCLUSIONS 

The vortex eigenfrequencies and magnetization profiles 
were calculated on the basis of analytical and numerical 
approaches. The low- lying part of the spectrum of spin 
excitation over the vortex ground state consists of dis- 
crete eigenfrequencies of magnetostatic origin. The cor- 
responding dynamic eigenmodes are localized near the 
dot center (the lowest translation mode) or are excited 
mainly outside the vortex core (radially symmetric mag- 
netostatic modes) . The model of the shifted vortex with 
no side surface charges^ well explains the results of our 
micromagnetic numerical calculations of the translation 
mode eigenfrequency. The set of magnetostatic eigen- 
functions Hn{p) and eigenfrequencies w„ was calculated 
as a function of the dot aspect ratio L/R. The first mag- 
netostatic eigenfunction can be approximately described 
as a uniform mode (no nodes). Whereas the high fre- 
quency eigenfunctions have a number of the nodes pro- 
portional to their number. The frequencies of the radial 
eigenmodes are proportional to y/L/R for thin magnetic 
dots with dot thickness L rj L cx and dot radii of about 
1 /zm. 
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